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We show at the example of the matrix element between pion states of a twist-2, non-singlet 
operator that Wilson twisted mass fermions allow to compute this phenomenologically relevant 
quantitiy at small pseudo scalar masses of 0(270 MeV). In the quenched approximation, we 
investigate the scaling behaviour of this observable that is derived from a 3-point function by 
applying two definitions of the critical mass and find a scaling compatible with the expected 
0(a 2 ) behaviour in both cases. A combined continuum extrapolations allows to obtain reliable 
results at small pion masses, which previously could not be explored by lattice QCD simulations. 
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1. Introduction 

Precision experiments allow to determine the parton distribution functions f(x, Qq) as a func- 
tion of the momentum fraction x and the energy scale <2q. At a fixed value of <2q, this function is 
usually parametrized as f(x) «i a (l — x)$ with exponents a, /3 that are obtained from global fits 
of the world-wide available experimental data (see ref. [[T]] for a general introduction). The bridge 
to results from non-perturbative calculations in lattice QCD are the moments J dxx n ~ l f{x) of the 
parton distribution functions which are directly computable with lattice technues. From the phe- 
nomenological side, there has been a lot of work to estimate the errors for these moments [3], 
These investigations also demonstrated clearly that the uncertainties of the parton distribution func- 
tions can easily reach 100%, in particular for the gluon distribution function. Hence, results from 
lattice calculations are highly welcome for phenomenologists and experimentalists. 

Lattice calculations for moments of parton distribution functions have progressed very much 
over the last years. We control the continuum limit [Q] and the finite size effects |^j, we understand 
fully the non-perturbative, scale dependent renormalization ^ of the operators needed to compute 
the moments and we can obtain matrix elements of such operators with a good precision. Nev- 
ertheless, even in the quenched approximation, there remains an important and hard to quantify 
systematic error: so far simulation data are obtained at only rather large values of the pseudo scalar 
mass of 600MeV of even higher. Clearly, an extrapolation to the physical value, i.e. a pion mass 
of about 140MeV, is difficult since contact to chiral perturbation theory cannot be established with 
such large pseudo scalar masses. 

In fig. 9 of ref. [JsJ] we give an example of the situation of our earlier calculation of a twist-2, 
non-singlet operator in pion states, correponding to the average momentum of the parton. The gap 
between the smallest pseudo scalar mass of about 600MeV reached in this work and the physical 
point at 140MeV appears to be quite large, leaving a substantial uncertainty in confronting the 
lattice calculations with the experimental results. 

In this contribution we will report about a way to overcome the uncertainty of the chiral ex- 
trapolation by employing Wilson twisted mass lattice fermions [Top which allows to regulate 
unphysical, small eigenvalues of the Wilson-Dirac operator, and hence can reach pion masses of 
well below 300MeV and thus provide a tool to use chiral perturbation theory to perform the last, 
small step to extrapolate to the physical point. See [11, [12|] for recent reviews. 



We emphasize that it is the main goal of this contribution to show that the gap in the pseudo 
scalar masses, seen in fig. 9 of ref. [§J, can indeed be closed when using Wilson twisted mass 
fermions. It is not the aim of the paper to provide a phenomenological number of the matrix element 
under consideration here and to compare to experiment. At the example of a 3-point function, the 
paper shall provide a demonstration that the Wilson twisted mass approach can be used to simulate 
small pion masses of 0(270 MeV) while keeping 0(a 2 ) cut-off effects under control as has been 



shown already for quantities that derive from 2-point correlation functions Q13| , 14) , |15Q . A more 



detailed account of this work will be published in ref. [|1 6|]. 
2. Lattice action and operators 



In this paper we will work on a lattice L 3 xT with twisted mass fermions the action of which 
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can be written as 

S[U,Y,W\ =a A Y J V{ x ){ D w + m + ipLY5'C?,)v{x) , (2.1) 

X 

with the standard Wilson-Dirac operator D\y and niQ, [i the bare un-twisted and twisted mass pa- 
rameters, respectively. Here and in the following y( x ) indicates a flavour doublet of quarks. 

In the present work, we will employ two definitions for the critical mass with the final aim of 
a combined continuum extrapolation of the results obtained in the two cases. The first definition 
of the critical mass is the point where the pseudoscalar meson mass vanishes, the second, where 
the PCAC quark mass vanishes. In the following we will refer to the first situation as the "pion 
definition" and to the second situation as the "PCAC definition" of the critical point, as introduced 
already in [15]. Both definitions should lead to f?(a)-improvement JT0| , 17, 18, 19], but they can 



induce very different 0(a 2 ) effects, in particular at small pseudoscalar meson masses [pCQ. 

The towers of twist-2 operators related to the unpolarized structure functions have the follow- 
ing expressions 

OjU^to = 2^T¥W7{ M1 D^} ^ a W(x) , (2-2) 

where {• • • } means symmetrization on the Lorentz indices, D^=D^ — D^i and = ^[V^ + 

V* ] . The flavour structure is specified by the Pauli matrices z a where we include here also the 
identity with T 4 = 2 • 1. In general one should perform an axial rotation in order to have the 
expressions for the twist-2 operators for the twisted mass formulation. However, the operators in 



eq. (2.2) with flavour index a = 1,4 do not rotate. We concentrate in this work on the twist-2 quark 
operator related to the lowest moment of the valence quark parton distribution function in a pion. 
In particular we will use for the up quark (the down quark can be treated in the same way) 



1 



1 



3 



(1 + T 



3> 



G\\ to = ~ W) 74 Da - t £ 7k D k y(x) . (2.3) 

J k=i 

The matrix elements of this operator can be computed in the standard way, described in refs. 
ill [22|]. We indicate with P ± (x) = y(x)y 5 i ^y{ x ) an d ^ = T '^' T " the interpolating operator for 
the charged pion. The ratio of the 2-point and 3-point functions 

C P (x 4 )=a 3J £(P + (0)p-(x, M )) , C 44 (y 4 )=a 6J £(P + (O)0 44 (y,y A )p-(x,T/2)) , (2.4) 

is related to the matrix element we are interested in. In particular if we perform a transfer matrix 
decomposition and we define R(y 4 ) = c^jji) > we obtain, in the limit when only the fundamental 
state dominates (0 <C y 4 T/2), the relevant bare quantity as 

Mbare = J_ -jR _ (2 . 5 ) 

m n 

The matrix element obtained in this way has to be renormalized by a multiplicative renormal- 
ization factor. To this end, we took the Zrgi as computed in ref. [|| 0] using the Wilson action. 
The renormalized matrix element has a well defined continuum limit and in the phenomenological 
relevant MS scheme is given by 

(x) MS (n,r m n ) = lim K ' RGU > pL=2GeV, (2.6) 

where Zrgi and/ MS (jU) were computed in ref. [S] (see this reference for further details). 
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3. Numerical results 



Our quenched simulations were performed for a number of bare quark masses in a corre- 
sponding pion mass range of 270 MeV < m K < 1 .2 GeV using the Wilson plaquette gauge action, 



employing periodic boundary conditions for all fields. In ref. [16] detailed tables with the numeri- 
cal results can be found. We performed simulations at fixed values of the twisted mass parameter 
jU while keeping mo to be the critical quark mass as obtained from the pion or PCAC definitions. 

Using the force parameter r$ = 0.5 fm to set the scale we determined the physical value of the 
pion mass and interpolated the values of (x) for the chosen values of rapsro. These values are close 
to the simulated ones, and so even a linear interpolation is usually sufficient. The lowest mass that 
can be reached corresponds to mps = 272 MeV. Note that we corrected our values for (jc) for finite 
size effects using the results of ref. 

In fig. |T] we show the continuum extrapolation of (jc), already converted to the MS scheme 
as explained in the previous section, at fixed pion mass for a wide range of pion masses. The 
scaling of (x) is in agreement with the expected 0(a 2 ) cut-off effects for both definitions of the 
critical mass. The slope of (x) as a function of a 2 appears to be rather small, at least for the pion 
masses investigated here. It also seems that the slope changes sign from positive to negative values 
when the pion mass is decreased. The plot also reveals that the continuum extrapolations for (x) 
performed independently for both definitions of the critical mass nicely agree. 




. m =368 MeV 
, 1^=368 MeV 
, m =728 MeV 
, mp S =728 MeV 
m =1051 MeV 
. 111^=1051 MeV 
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Figure 1: Left panel: (x) MS as a function of a 2 for different values of the pion mass. Right panel: (x) MS (/z = 
2 GeV) extrapolated to the continuum as a function of the pion mass. Open squares represent results that are 
obtained from a combined continuum extrapolation of earlier Wilson and clover- Wilson simulations, ref. M\, 
The filled cicles represent our results using Wilson twisted mass fermions. The open circle denotes a result 
which is not corrected for finite size effects and the diamond corresponds to the experimental point. 

In fig. |l| (right panel) we show (x) MS in the continuum as a function of the continuum pion 
mass. The open squares represent our values obtained earlier from a combined continuum extrap- 
olation of Wilson and clover-improved Wilson fermions using the Schrodinger functional scheme 
[Bl. As usual, such simulations have to stop at a pion mass of about 600 MeV. At such high pion 
masses it becomes very difficult, if not impossible to compare the simulation results, even when 
extrapolated to the continuum limit as done here, to chiral perturbation theory. 
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The right panel of Fig. [I] shows that with Wilson twisted mass fermions, represented by the 
filled circles, the large gap between a pion mass of about 600 MeV, as the lower bound for standard 
simulations, and the physical value can be bridged. The Wilson twisted mass results for (x) are 
consistent with a linear behaviour. Note that we add in the figure a point denoted with an open 
circle which is not corrected for finite size effects and is shown for illustration purpose only that 
also pion mass below 300MeV can be reached with Wilson twisted mass fermions. 



4. Conclusions 



In this contribution we have computed the lowest moment of a pion parton distribution fun- 
tion. In particular, we concentrated on a twist-2, non-singlet operator the matrix element of which 
corresponds to the average momentum of a parton in the pion. We demonstrated, that also for 
such matrix elements pseudo scalar masses below 300MeV can be reached. At the same time, 
the cut-off effects even at such small pseudo scalar masses appear only in 0(a 2 ) and remain small 
and controllable. This clearly opens the door to make contact to chiral perturbation theory and to 
perform reliable extrapolations to the physical point with a pseudo scalar mass that corresponds to 
the experimentally observed pion mass. 

In the present work we have employed two definitions of the critical mass, the pion and PCAC 
mass definition [14, 15]. We could perform a controlled continuum extrapolation of (x) at least 
down to pion masses of about 270 MeV by combining the data obtained with the two definitions. 
Of course, in principle, also overlap simulations are able to reach such values of the pion mass. This 
will come, however, at a much higher simulation cost [23]. The right panel of our final figure, fig. [j], 
clearly demonstrates that with our present setup it is possible to go beyond standard calculations 
and bridge the gap between large pion mass values of 600 MeV and the physical value of the pion 
mass. 
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